Training-induced criticality in martensites 
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We propose an explanation for the self-organization towards criticality observed in martensites 
during the cyclic process known as 'training'. The scale- free behavior originates from the interplay 
between the reversible phase transformation and the concurrent activity of lattice defects. The 
basis of the model is a continuous dynamical system on a rugged energy landscape, which in the 
quasi-static limit reduces to a sandpile automaton. We reproduce all the principal observations 
in thermally driven martensites, including power-law statistics, hysteresis shakedown, asymmetric 
signal shapes, and correlated disorder. 
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Experiments in martensites reveal intermittent behav- 
ior with power-law statistics [ij, |2, Isl, \^, showing an 
intrinsic complexity comparable to that of turbulence, 
earthquakes, internet networks, and financial markets. 
Criticality is known to be an issue of great significance 
in contemporary science, giving a framework for under- 
standing the emergence of complexity in a variety of nat- 
ural systems [5|, Ig . Within materials science, criticality 
has been recognized in the last years as a key factor in 
crystal plasticity, brittle fracture and damage [7| . In the 
present paper we develop a model explaining why similar 
behavior is also observed in martensitic transformations. 

Reversible martensitic transformations involve a coor- 
dinated distortion of the crystal lattice and belong to the 
class of ferroelastic first-order phase changes with ather- 
mal character @, 0. In such systems the macroscopic 
strain discontinuity typically splits into a set of bursts 
(avalanches) corresponding to transitions between neigh- 
boring metastable states. The individual avalanches can 
be detected through the measurement of the intermit- 
tent acoustic or calorimetric signals. The size distribu- 
tion observed in shape-memory alloys (Cu-Al-Ni, Cu-Zn- 
Al, Cu-Al-Mn, Ni-Mn-Ga) was shown to be scale free 
[l|, y, y, y, y . Despite the apparent similarity with driven 
ferromagnetic systems, where the scale-free Barkhausen 
noise has been known for a long time, the experiments 
on memory alloys show features not observed in magnets, 
and which are instead reminiscent of plastic shakedown. 
In particular, the critical character of the avalanches [3,|3[ 
and the smoothing of the hysteresis profile [3|, llO| emerge 
only after multiple thermal cycling through the transi- 
tion. 

The mechanism leading to training-induced critical be- 
havior in martensites strongly resembles the phenomenol- 
ogy associated with self-organized criticality (SOC) [5|. 
The SOC paradigm in the form of a sandpile automa- 
ton has been applied to martensitic transformations in 



exploited the similarity between martensites and mag- 
netics by interpreting both in an Ising-type framework, 
with zero temperature and quenched disorder 1^ . In this 
context the power-law over a few decades of avalanche 
sizes is viewed as a sign of proximity of the system to 
a classical critical point. Criticality then emerges only 
as a result of the external tuning of the disorder. Fur- 
thermore, the symmetric avalanche shapes with scaling 
collapse which are expected in these models [1^ , contrast 
those experimentally recorded [13[. More recent model- 
ing has focused on the direct simulation of martensitic 
transformations within the framework of elasticity the- 
ory |l4l |. While the corresponding numerical tests show 
some scaling in avalanche sizes, the system is unable to 
memorize its state of disorder upon unloading, failing 
therefore to exhibit the effects of training. 

A key experimental observation left aside in the pre- 
ceding theoretical work is the dislocational activity assist- 
ing the development of the phase transformation in these 
materials [lO|, ll5|, UM, llll, ll8| . It has been repeatedly ob- 



lll : which however, lacked a connection to the physics 



of martensitic transformations. A different set of models 



served that dislocations are indeed introduced in shape 
memory alloys during training. In particular, under pe- 
riodic driving the degree of defectiveness first increases 
monotonically and then saturates [16|, ll7|. Our model 
shows that this dislocational activity is highly correlated 
and is ultimately responsible for the scale-free character 
of the reversible behavior of martensites. More precisely, 
the attainment of criticality is due to the ability of the 
crystal to develop an optimal amount of disorder. 

To describe in an unified way the processes involv- 
ing both dislocations and phase boundaries, we consider 
a prototypical 2D system of kinematically compatible 
elastic units resulting from a suitable triangulation of 
a square lattice (see an example of such procedure in 
llOf ). To each unit with index i we assign a multi-well 
strain energy function depending on a single scalar or- 
der parameter e which, in turn, is a combination of the 
components of the discrete strain tensor. The adiabatic 
elimination of the harmonic non-order-parameter strain 
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FIG. 1: Schematic representation of the periodic potential 
/(e; r) with the shaded box highhghting a single domain of 
periodicity. Insets on top show the lattice structures corre- 
sponding to the bottoms of the potential wells, with lattice 
cells marked for the austenite A and for the martensite vari- 
ants M±. The austenite A' corresponds to a different period 
of the potential. 



variables by means of the equilibrium equations and the 
kinematic compatibility constraints leads to a non-local 
elastic energy of the form 



$(e;T) = ^/(e,;T) + i^X, 
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where K = {Kij} is the kernel of the long-range elastic 
interactions, and / is a periodic function as in Fig. [T] 
In each period we use the three-parabolic approximation 
f{e;T) = i(e - w)^ + g(T)s^, where s = in the high 
symmetry phase (austenite) and s — ±1 in the two vari- 
ants of the low symmetry phase (martensite). The pa- 
rameter w = d + es, where e is the transformation strain, 
defines the location of the bottoms of the energy wells; 
the integer-valued parameter d specifies the period of /. 
We emphasize that no randomness have been assumed in 
the model. 

The global periodicity of the energy takes large shear- 
ing distortions into account [20[, so that both the phase 
change and dislocation formation can be handled simul- 
taneously. When the transformation strain is small and 
the energy barriers for slip are much higher than the bar- 
riers for the phase transition and twinning ('weak tran- 
sitions'), no lattice-invariant shears occur (as with the 
R-phase of NiTi [10| , or the pre-martensitic transforma- 
tion in Ni-Mn-Ga [J|) and modeling can proceed accord- 
ing to the Landau theory. Such phase changes do not 
generate significant dislocational activity and are largely 
reversible, as is assumed, for instance, in lj,l2l|. On the 
contrary, in 'reconstructive transformations' the trans- 
formation strain is large, as it lays at the boundary of 
the periodicity domain (e.g. the ideal Bain transforma- 
tion from bcc to fee). In those cases the energy barriers 
to slip are only as high as the transformation barriers, 
and the phase change advances 'half way' towards the 



formation of dislocations [l9[ . As a result, defects prolif- 
erate making the transition irreversible |22l | . In between 
these two extremes a range of possibilities exists, where 
defect formation plays an increasing role as the transfor- 
mation strain gets closer to the boundary of the maxi- 
mal periodicity domain. In particular, all the martensitic 
transformations considered in the experiments that we 
are concerned with [ij, [3, y, Ij, Isf , involve a transforma- 
tion strain e which is very close to the ideal Bain strain 
[23|. When slowly driven, these systems are expected to 
exhibit cell deformations not confined to one periodicity 
domain, but rather extending to an unbounded portion 
of the periodic energy landscape. This involves forma- 
tion of dislocations, which our kinematic compatibility 
assumption does not exclude, as is exemplified in Fig. [21 
where we compare the standard representation of a dis- 
location (b) with the one adopted in this paper (a). 





FIG. 2: Two representations of a dislocation, (a) Continuous 
lattice deformation involving partial slip, as is interpreted in 
the text, (b) The same atomic configuration viewed within 
the classical interpretation involving a discontinuous deforma- 
tion and a non-zero Burgers vector. The shaded units refer 
to austenite in the wells A and A' in Fig. [T] 



We drive the system quasi-statically trough the func- 
tion g{T), by changing the temperature T. Since the 
transformations are typically athermal [9|, we consider 
a purely mechanical setting with overdamped dynamics. 
By making this assumption we treat the acoustic emis- 
sion as a part of dissipation. In the limit of infinitely slow 
driving, the dynamics projects on the local minima of 
the total energy $, which form a discrete set of branches 
e = s{g), with e~ < e < e'^ , where the extremes e 



correspond to marginally stable configurations 2J|. For 
piece- wise parabolic / as in Fig. [U the limits e^ of each 
branch can be explicitly written as a function of g and 
e. When such limits are reached, the instability resolves 
through a fast event (avalanche) which brings the system 
to another equilibrium branch. In this way the dynamics 
becomes piece- wise continuous (see l25| for a study of the 
ID case). 

We proceed by eliminating through minimization the 
linear elastic strain e at given d and s obtaining e = 
(1 + K)~^w. The relaxed energy is of the Ising type 
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where the two discrete spin variables describe the phase 
transformation (s) and the plastic slip (d), and J = 
(1+K)^^ — 1. Since the energy ^ is supposed to penalize 
the inhomogeneity of the field Wi induced by either phase 
boundaries or dislocations, we assume that the corre- 
sponding term has the particular form | J2 Jij (^i ~ ^i)^ 
so that Jii = ~J2j=iiJij- We furthermore assume the 
kernel J to be of the ANNNI type, to account for the 
competing interactions driving both the coarsening and 
the refinement of the microstructure 2l|, |26|. Specifi- 
cally, we consider Jy = Ji > for nearest neighbors, 
and Jij = — J2 < for ncxt-to-ncarcst neighbors; the 
ensuing diagonal-dominated structure of the matrix J 
does not prevent the original matrix K in ([1]) from be- 
ing dense. Suitable inequalities ensure the metastability 
of the individual equilibrium branches [25|; in particular, 
the condition Ji > 2 J2 is sufficient for our automaton to 
reach a steady state (we use Ji = 0.062, J2 = 0.03 in the 
simulations presented below). 

Under the above hypotheses, the piece- wise continuous 
dynamics becomes a sandpile automaton, whose main 
variable is the elastic strain 6i = Ci ~ Wi representing 
the 'local height'. Once a cell becomes unstable (the 
condition e~ < e^ < e+ is violated), 5i is updated as 
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S, - 4(Ji - J2)r 
J ' Sj -\- JiT, j nearest neighbors of i. 



(3) 



Sk ^ J2T, k next-to-nearest neighbors of i, 



where r = ±e for phase transitions, and r ~ ±(1 — 2e) for 
slips. Since each update may make new sites unstable, 
the updates continue at constant g until the system is 
fully equilibrated. For J2 = 0, e = ^. and 5 = the 
automaton ([3]) reduces to the standard BTW sandpile 
i. 

We implement the model numerically on a 501 x 501 
grid for an almost reconstructive transformation with 
e = 0.47 and open boundary conditions. The initially 
homogeneous austenite is chosen to contain only a min- 
imal dislocation loop. The crystal is then thermally cy- 
cled through the complete transformation, g being a pe- 
riodic triangular function of computational time. The 
intensity of the acoustic bursts registered experimentally 
is linked to the size of the avalanches (total number of 
updates before stabilization). Fig. [3^ shows the devel- 
opment of the phase microstructure during the training 
period. The level of plastic deformation is monitored 
through the density p of nearest neighbors with differing 
values of di. Fig. [8)3 shows the formation of dislocations 
induced by training and marked by the steep initial in- 
crease in the variable p (Fig. [4^). The creation of cor- 
related dislocation microstructure quickly saturates, in 
accordance with the experiments [l6|, ll7[. In Fig. |4)d, we 
observe the smoothing effect of the self-organized defects 
on the cooling curves (and hence on the hysteresis cy- 
cle). The dislocational activity leads to the increase of 




FIG. 3: Evolution of the phase and defect microstructures in 
the lattice during thermal cycling. A minimal dislocation loop 
was placed at the middle of the system, (a) Configuration 
of the martensitic phase domains, represented by the field 
Si, after cycle 1 and cycle 2000 (turquoise and blue indicate 
s = 1 and s = — 1, respectively), (b) The corresponding 
configurations of the slip variable d. White and brown colors 
indicate d = G and d 7^ (mostly |d| = 1), respectively. The 
black contours separate elements with different d. 



the martensite starting temperature, similarly to what 
is reported experimentally [3| . The parallel development 
of criticality is indicated by the emergence of the power- 
law statistics for the avalanche sizes (Fig. [3J:). At the 
beginning of the training period the avalanche distribu- 
tion is supercritical with a peak at large sizes evident 
from the sharp initial cooling curves in Fig. [4)d. The 
eak eventually vanishes (Fig. lit), as in the experiments 
Two further predictions of the model matching ex- 
perimental data [l3| concern the strong asymmetry of 
the avalanche shapes, and the absence of their scaling 
collapse (see Fig. [4jl). Similar effects are also observed 
in Barkhausen noise, earthquakes, and dynamic fracture 
(see the discussion in|27|). Fig. |4h, shows the distribu- 
tion of avalanche durations (number of simultaneous up- 
dates in an avalanche) predicted by the model, which 
deviates from a power-law. This indicates that in the 
present framework a scale-free size distribution does not 
always imply scaling in time. While this prediction for an 
idealized system (with neither inertia nor fixed pinning 
sites) disagrees with the scaling for durations reported 
experimentally [l|, our model does generate an almost 
power-law distribution of durations after the introduc- 
tion of a small amount of quenched disorder represented 
by a Gaussian distribution of 5i with zero average in the 
initial configuration. Such modification, however, does 
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FIG. 4: Evolution to criticality and features of the critical 
state, (a) The dislocation density p during the first 2000 cy- 
cles, (b) Cooling curves representing phase fraction vs. tem- 
perature in arbitrary units (cycle 1, solid line; cycle 200, 
dashed line). The corresponding experimental data [3|] are 
shown in the inset for comparison (1st and 24th cycles), (c) 
Steady-state power-law distribution of the avalanche sizes, 
and non-power-law distribution of durations (displaced one 
decade lower for clarity), (d) Asymmetric avalanche shapes. 
The inset shows the corresponding experimental data [13| . 



not influence the power-law distribution of avalanches, 
nor does it affect the pulse asymmetry. 

When in the original setting with two variables s and 
d the phase transformation is suppressed (no variable s), 
the model describes the micromechanics of stress-driven 
intermittent plastic flow in crystals [28|. In this case the 
system is defined by the single integer-valued order pa- 
rameter d as in the phase-field description of plasticity 
[29[, to which our model offers an analytically accessible 
automaton alternative 



32l |. What is more, the present 



scheme of plasticity does not require ad-hoc procedures 
for the nucleation and annihilation of dislocations, as in 
discrete dislocation dynamics [30[. When, in contrast, 
the variable d can be neglected in the original setting, 
as in weak martensitic transformations (or in magnet- 
ics), the model still generates several decades of power- 
law avalanche distribution in accordance with the exper- 
iments in Ni-Mn-Ga [J| , but only in the presence of some 
quenched disorder represented again by a Gaussian dis- 
tribution of Si in the initial configuration. In this case 
the emergence of limited scaling does not involve train- 
ing, and can be explained b y th e proximity of the system 
to a classical critical point [12l |. 

In summary, the proposed model accounts for all the 
main observed phenomena accompanying the training 
process in martensites leading to criticality. The agree- 
ment with experiment clearly indicates that SOC origi- 



nates in these systems as a result of the interplay between 
the reversible phase change and the irreversible develop- 
ment of an optimal amount of plastic deformation. 
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